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A numerical model of the Ares I upper stage main propulsion system is formulated based 
on first principles. Equations are written as non-linear ordinary differential equations. The 
GASP fortran code is used to compute thermophysical properties of the working fluids. 
Complicated algebraic constraints are numerically solved. The model is implemented in 
Simulink and provides a rudimentary simulation of the time history of important pressures 
and temperatures during re-pressurization, boost and upper stage firing. The model is 
validated against an existing reliable code, and typical results are shown. 
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temperature 
heat- transfer area 
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Superscript 
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viscosity 

density 


ambient 
engine 
film layer 
pressure gas 
inlet 

propellant liquid 

interaction through common bulkhead 
wetted tank surface 
tank or total 

propellant or oxidizer vapor 
tank wall 


derivative with respect to time 


I. Introduction 

T he fault detection, diagnosis and recovery laboratory (FDDR) at Marshall Space Flight Center will 
provide simulation of key subsystems of the Ares I launch vehicle and test innovative fault detection 
algorithms for possible employment on the vehicle: Since this development is proceeding in parallel with 
the actual vehicle design and construction, little or no actual test data exists. That is, actual time histories 
of pressures, temperatures, accelerations and other measurements do not yet exist. In order to provide a 
realistic simulation of the vehicle and environment, models based on physics must be developed. This work 
focuses on the Ares I upper stage main propulsion system (MPS). 

In the past 50 years, many researchers have sought numerical models of pressurized propellant and 
oxidizer tanks. These models attempt to capture the heat and mass transfer behavior of the tank and 
working substance which for cryogenic fuels is multiphase due to boil off. As early as 1965, Roudebush 
obtained numerical solutions accounting for axial temperature distributions in the ullage gas and tank wall. 1 
Pasley 2 proposed an ODE model, with continuous ullage gas flow and demonstrated good comparison with 
hardware test data from the Surveyor spacecraft. Kendle 3 compared experimental results with models 
assuming both mixed and stratified ullage volumes. Bush and associates 4 derived models for both ‘small’ 
and ‘large’ ullage cases of a closed tank. For the small ullage case, the time history equations were reduced 
to a non-linear algebraic set. For the large ullage case, phase change complicated the problem and resulted 
in a set of 1-D PDEs. 

More recently, several researchers have focused on modeling the X-34 MPS. 5-9 X-34 uses RP-1 (kerosene) 
fuel and liquid Oxygen (LOX) oxidizer. Mcdonald et al 5 performed CFD flow analysis of the X-34 MPS 
to detail fuel and oxidizer eg time history, tank bottom pressure history, and fuel/ullage stratification, a 
paramount constraint since the X-34 is launched in a horzontal attitude. Their results provide qualitative 
insight into the current solution. Majumdar and Steadman 6 propose a finite volume numerical model in- 
cluding second law effects. They have developed a software tool called Generalized Fluid System Simulation 
Program (GFSSP) for network flow analysis. GFSSP has been validated against the Propulsion Test Arti- 
cle Helium Pressurization System 7 at NASA Stennis, and the X-34 MPS. 8 Thus it is considered a reliable 
simulation for the problem at hand and will be used to validate the models developed herein. 

Here we capitalize on a model due to Estey et al. 12 We are not concerned with spatial effects, hence the 
model is formulated by choosing a few control volumes, and treating each as a lumped element. In section 
II, we detail the system equations and correlations. Section III describes specifics about the Simulink imple- 
mentation. Section IV presents the results of validation and nominal performance for boost and mainstage 
operation. 
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II. Development of the State Space Model 

The Ares I upper stage fuel and oxidizer tanks are depicted in Figure 1. Note that all relevant heat 
and mass transfer relationships are depicted, however, unique subscripts are not assigned for each tank, so 
the terms must be interpreted in their context (within each tank). The current design employs a common 
bulkhead between the Hydrogen and Oxygen tanks, and submerses the cryogenic Helium tanks in the H 2 
tank, thus requiring additional heat transfer interactions which are depicted. Each tank is divided into 
four control volumes, namely, the tank wall, ullage gas volume, film layer, and fuel/oxydizer liquid region. 
Temperature dynamics of the common bulkhead are not modeled-the bulkhead is simply treated as a thermal 
resistance. 



Figure 1. Arrangement of Ares I Upper Stage fuel and oxidizer tanks with heat and mass transfer interactions depicted. 


II. A. Liquid Oxygen Tank 

The Liquid Oxygen Tank is modeled first since it is pressurized in all modes with gaseous Helium. Estey 12 
provides a state space model for a single substance pressurized tank in blowdown mode. The Ares I MPS 
will use Helium pressurant in both liquid Oxygen and liquid Hydrogen tanks. Thus we must expand Estey’s 
model to accomodate a gas/vapor mixture in the ullage space, and addition of ullage mass during tank 
re-pressurization and fuel expulsion. As shown in Figure 1, each tank is divided into four control volumes. 
The conservation of energy is then written for each control volume as follows. For the ullage space, the 
conservation of energy equation is 

Qg “1“ Qv QH2 "P fflyRyTy ~b ^Xg RgTg {iXXyCy T TTLqCg ) Ty (l) 

The film layer is considered infinitessimally thin, and thus^ we obtain the same equation as Estey 12 

Ql ~ Qv + m v hi ~ rn v h v = 0 (2) 

The term m v (h v — hi) can be rewritten as 

rh v [(h v - h vs ) - (hi - h is ) + (Ks - his)] = m v [{c v + R V ){T V - T f ) - ci{Ti - T f ) + h fg ] (3) 
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(4) 


For the liquid region, the first law equation is unchanged from Estey 12 

Qp-Qi + W = m e {v 2 e /2) + mciTi 

The tank wall conservation of energy equation is also unchanged from Estey 12 

Qw ~ Qp Q g — TTI"uj C yj ( 5 ) 
The heat transfer terms above are expanded as 

Qg = hgA g (T w —T v ) ( 6 ) 

Q v = hfAf{T f — T v ) ( 7 ) 

Qi = h s A f (Ti — Tf) ( 8 ) 

Q p = h t A s (T w —Ti ) ( 9 ) 

Qw = h w A w (T a -T w ) ■ ( 10 ) 

Qo 2 — ho 2 Ao 2 {T v -T h 2 )- ( 11 ) 

The heat transfer coefficients h g and h s = hf are determined from the natural convection correlation 

hi = K H C(k f /l s )X n 

where 

X - (Gr)(Pr) 

Gr = (iy f gl3,\&T\/fi) 

Pr = c p ffif/kf 

Gr and Pr are the Grashof number and Prandtl number of the ullage gas, respectively. According to Ring, 15 
C = 0.27, n = 0.25, and Kh = 1.. The length scale l s is set to the diameter of the tank. The volumetric 
thermal expansion coefficient is defined as 


Pf = ~ 


1 

Pf 



and must typically be derived from finite differencing of p wrt T. All thermophysical and transport properties 
required to compute Gr and Pr are found from ullage pressure and temperature using a Matlab executable 
(mex) of the Fortran based GASP code. 16 In the O 2 tank, the ullage gas is assumed to be a mixture of 
Oxygen and Helium, and the properties are computed by a weighted average of the pure substance properties 
based on mass ratios. 

Substituting the heat transfer terms and Eqn 3 into Eqn 2, we obtain 

m v [(cv + Rv)(Tp — Tf) — ci(Ti - Tf) + h/ g ] = h s Af(Ti - Tf) + hfAf(Tf — T v ) (12) 


which can be solved independently for the rate of evaporation rh v . Substituting the heat transfer terms into 
Eqn 5, we get 

m w c w T w = h w A w (T a - T w ) - h t A s (T w - T L ) - h g A g (T w - T v ) (13) 

which can be solved separately for the time rate of change of the wall temperature. The remaining equations 
are coupled during fuel expulsion due to the boundary work exerted by the ullage gas mixture on the liquid. 


II.A.l. Pressurization ( Fuel Expulsion) Model 

Assuming a common temperature in the ullage space, and applying Dalton’s law, the total ullage volume is 


V to t = 


m v R v m q R 


Pv 


+ 


Pa 


T„. 


( 14 ) 
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The boundary work modeled as quasistatic work of expansion, W = P t (dV g /dt ), where P t is the total ullage 
pressure. Differentiating V tot yields 


dVtot ( TTly R v TTlyRsyPy t TTTgRg TYlgRgPg'^ ^ | f TYlyRy | TTlgRg^ ^ 

-pT + -p~. - —pT ) v + ~pPJ v ' 


Then 


W = (P v + P g ) 


( TTlyR^y TYlgRr 

[~pT + 


p„ 


Ty + 


TTlyRy TYlyRy Py TTlgRg TYlgRgPg 


P 2 + Po 


P 2 

9 


Ty 


1 9 J \ Py 1 v 1 g 

Substituting for the work term and the heat transfer terms in Eqn 1, we obtain 

{ YflyCy + m g Cg)fy — ( TTlyRy + TYlgRg)Ty HgAg{T w Ty ) — h0 2 g (Ty — T fj 2 ) + H f A f (T f — Ty ) 

- (Pv + Pg) 


/TTlyRy TYlgRg\ • I TTly Ry TTlyRyPy | TTlgRg TYlgRgPg ^ ^ 

w~ + ~p^~~^r) v 


(15) 


(16) 


(17) 


Likewise, Eqn 4 becomes 


-m e vl + miciti = h t A s (T w - T t ) - h s A f (Ti - T f ) - ho 2 iAf(Ti - Th 2 ) 


+ (Py + Pg) 


ITTyRy ^ TTlgRg 


Py 


Pn 


T v + 


TYlyRy TTlyRyPy ^ TTlgRg TYlgRgPg 

~Ty W ' 


+ 


Pn 


PI 


T. 


The tank total volume is invariant, thus given 

TlTtyRyRy TTlg RgTy 


Vt=* 


Py 


■ + 


Pn 


+ mivi 


and differentiating, we obtain 


Q TTly n e l Ry Ty ^ TflyRyTy TTly RyTy Py ^ TYlg ne ^RgTg ^ TTlg RgTy TTlg RgTy Pg 


Pn 


Py 


R 2 


Pn 


Pn 




(18) 


(19) 


( 20 ) 


Where Th vn et and Tn gnet include the sum of evaporation, venting and addition of ullage Oxygen vapor and 
Helium gas respectively. The system of equations is completed by enforcing an isothermal constraint in the 
ullage space. 

T 9 = ^= T v = ^- (21) 


TTlg Rg 


TTly Ry 


Solving for P g and differentiating, we obtain 


p _ p m g R g | m g Rg m g rhyRg 
9 V TTlyRy TTlyRy TTlyRy 

For ease of algebraic, or numeric solution, Eqns 17, 18, 20 and 22 are arranged into the matrix form 


(22) 


" -Gil 

G12 

0 

— G14 


' Pv ' 


' # 1 1 

— G21 

G22 

0 

-G24 

< 

Ty 

— 

h 2 { 

G31 

-Gi 2 

G33 

G34 


Ti 


h 3 ( 

. G41 

0 

0 

-1 


{ Pg J 


l h 4 J 


(23) 


Where 


On = 


TTly RyTy 

P.? 


TTly Ry TTlg Rg 
O' 12 = — — H 


Pn 


Pn 


(24) 

(25) 
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G*14 

a? ST 

II 

(26) 

G 2 i 

— —m v R v T v 4- 

(27) 

g 22 

= myCt, + m g Cg + (Py + Pg)G 12 

(28) 

G 2 4 

II 

2 

to 

£3 

to 

+ 
ha 1 

50 tO| <2^ 

(29) 

G 31 

- (Py+Pg)G 11 

(30) 

G 33 

= mici 

(31) 

G 34 

= (Py + Pg)G 14 

(32) 

G 41 

m g R g 

(33) 


and 


rflgRgTg 

Hl ~ — p— 


Tfly RyTy 

Pv 


+ m e vi + Virhy 


\h(qp rp__ \ , u - A ./CP. T 1 \ PfjPv ...... PyPc, 


H 2 = h g Ag(T w - Ty) - ho 2 Ag(T v - Th 2 ) + hjAf(T f - T v ) - 


Pv 


-m v - 


r 9 


(34) 

(35) . 


H 3 = h t A s (T w -TO - h s Af(Ti - Tf) -hoAW -T H ,) + (P„ + p g )(^k + V^s) Tv - im et , e 2 (36) 

Hi= _n^ + rn s m^ 9 _ 
m v R v ml R v 

An algebraic solution of Eqn 23 is possible, but cumbersome, so we rely on a numeric solution for the state 
derivatives in Eqn 23. . 

II. A. 2. Re- Pressurization and Re- circulation Model 

During ground operations, and first stage boost, the liquid fuel mass inside the tank is essentially con- 
stant, and the boundary work term may be neglected. This modifies the previous equations in that terms 
(? 2 i) G 2 a, G 31 , G 12 , and G 34 are zero. Thus the equations for T v and T; may be solved directly as follows. 

T v = — - — [h g A g (T w — Ty) — ho 2 A l g(T v — Th 2 ) + k f A f (T f — T v ) + ( rhynetRy + m g netRg)Tv] 

' (38) 


rriyCy + m g Cg 
1 


Ti = [h t A s (T w - Tt) - h $ A f (Ti - 2» - ho 2 A?(T t -T m )~ - T u )] 


mici 


(39) 


Where mi and Tu are respectively the mass flow rate and return temperature of the recirculation loop. The 
equations for P g and P v remain coupled and can be written as the linear system 


1 — Qai 

G 14 G\\ 



Pg 1 

L_J 

< 

A J 

11 


Wig net Rg 771 <y lily net Rg 

TTXy Ry 777.^ Ry 


ZZhmxufc. 


RyT v _j_ m v R v T v _ _|_ Wgne.tRgTg rrigRgTy 

Ry Py Pg Pg 


(40) 


II. B. Liquid Hydrogen Tank 

The Liquid Hydrogen Tank is somewhat more complicated since its ullage gas is Helium for re-pressurization, 
and Hydrogen during engine firing. It is also more complicated due to heat transfer interaction with the 
cryogenic Helium tanks which are submersed in the liquid Hydrogen tank. 


II.B.l. Pressurization (Fuel Expulsion) Model 

During fuel expulsion, gaseous Hydrogen is added to maintain ullage pressure. The ullage gas will contain 
some He from re- pressurization, however we consider the He mass to be negligible in light of the rapidly 
growing Hydrogen vapor mass. The energy balance equation for the ullage volume is 

Qg T Qv QtJeg 'm v CyTy TYlyRyTy JllyU RyTy (4l) 
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The boundary work term for a single species ullage is 


W = 


Pv 


m v T v + mn v t v ) P v — m v T v P v ) R 


(42) 


where m v and rh u u are vapor mass added by evaporation, and ullage pressure regulation respectively. Thus 
substituting the heat transfer terms and boundary work, Eqn 3b becomes 


TRyCyTy {ply P tilull) RyPv hgAg(T W Ty} + hfAf(Tf Ty) hff e gA}J e g{T v T// e ) 

- ~ ( (rn v T v + m v t v ^j P v - m v T v P,)j R v (43) 

The film layer dynamics are identical to those of the Oxygen tank. The liquid Hydrogen control volume 
differs from that of Oxygen because of the single species work term, and the extra heat transfer term due to 
interaction with the cryogenic He tanks. 

ip 1 

Q P + Qh 2 ~ Qi — Qh&i A W = m e -j- + miCiTi (44) 

Which after the appropriate substitutions becomes 

vP 

1 ^ le ~2 m i c iPi = h t A s (T w — Ti) + hH 2 AH 2 {To 2 ~ Pi) ~ h s Af(Ti — Tf ) — hnei^-Hei{Pi — Tjje) 

+ - 5 - (( rh v T v + m v T v ^ P v - m v T v P v ) R v ( 45 ) 

Note that subscripts in these equations must be interpreted in their context, i.e. 7} immediately above refers 
to the temperature of liquid Hydrogen, wheras in subsection II. A., it refers to the temperature of liquid 
Oxygen. The tank wall energy equation is identical to Eqn 13. The total volume equation for a single 
species ullage is V t m m v R v T v / P v + rriivi which when differentiated yields 


Q AlyRyTy ^ TTly RyTy 


P« 


P, 


TTly RyTy P v . . Rx)Ty 

— m v vi- m e vi + m u u - 


P ,2 


Pv 


(46) 


To facilitate solution of the state derivatives, Eqns 43, 45 and 46 are arranged in the matrix form 


Tn>v Ry T?) 

p2 
1 V 

ni v P-v 
Pv 

0 


[ Pv \ 



m v R V T V 
Py 

— TTlyCy — rriyRy 

0 

< 

1 rp 1 

| lv 1 

> = < 

1 Hl 

TB.-ai.Tv T.v. 

L p v 

TTly Ry 

-mid _ 


l * J 

1 

l j 


where 


(47) 


Illy RyTy . RyTy 

Pi = 5 b TUyVi + m e vi - mull —p, — (48) 

V Ty 

H 2 —rhyTyRy ~ (ijly + lh U u) RyTy ~ hgAg(T W ~ Ty) ~ 1l f A f {T f ~ Ty ) + k H eg A^Ieg (Ty ~ Ttf e ) (49) 

Ip 

Hs = rh e -p- — rhyTyR v — h t A s (T w — Ti) — hn 2 AH 2 (To 2 — T{)+h s Af(Ti — Tf)+hHeiAnei(Ti — Tne){ 50) 


Clearly the LHS of Eqn 47 is a Hessenberg permutation, so the equations for P v and T v can be solved 
independently of T) . Also due to this structure, the equations can be algebraically solved relatively easily 
yielding 


Pv 1 

1 

(Cy + R v )Py 

Pv 


m v T v R v c v 

Wly TyCy 

Ty 

l 

Pv 

1 

r 

m v Cy 

m v c v 

T J 

1 

Py 

mi ci 

0 


0 

0 

1 

mici 



( 51 ) 
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II. C. Pressurant Subsystems 

Prior to launch, the cryogenic He bottles are pressurized to 3277 psi or 22.59 MPa, which is a supercritical 
pressure. The ten bottles have a total volume of 1.00 m 3 . The initial stored He mass equals 200 kg, which is 
based on the assumption that the cryogenic He initial temperature is equal to that of the liquid H 2 ( 20K). 
Prior to injection in the fuel and O 2 tanks, the He passes through a regulator which reduces its pressure to 
435 psi or 3.0 MPa. We simulate the pressurization control as a simple relay with on/off thresholds set to 
35 and 40 psi respectively for the H 2 tank and 25 and 30 psi for the O 2 tank. 

In order to determine the mass flow rate of ullage gas into the tank when the control valve is open, we 
use the relationships from Palm. 13 The critical pressure for an ideal gas is 

/ 9 \7/(7-i) 

p ‘ = Utt) p > < 52 > 


where Pi is the supply pressure. For typical gases, the critical pressure tends to be about 1/2 of the supply 
pressure or higher. Since the tank pressure is around 1/10 to 1/20 of the supply pressure in this case, we 
consider the flow into the tank to be sonic, and use the following relationship for mass flow rate 


'rii'uii — CdA 


Pi 

vm 


,7+1 


(7+l)/(7-l) 


(53) 


Where A is the orifice area, T\ is the supply temperature, and Cd is a discharge coefficient which is assumed 
to be unity. 

Knowing the net mass flow rate out of the cryo He tanks and into the fuel and O 2 tanks, we predict 
the time history of mass in the cryo He tank simply by integrating. The cryo He pressure history is then 
calculated by GASP 16 from current density assuming temperature equals 20K. 


II. D. Thermophysical Properties 

The film temperature is assumed to be the working fluid saturation temperature at current tank pressure. 
These are calculated from saturation pressure by numerical solution of the equation 

l°gio Psat — ji + ^ + hT + j^T 2 + jsT 3 + jaT 4 4- jjT 3 (54) 

where the j coefficients for Hydrogen and Oxygen are shown in Table 1. Since the equation above is a 
polynomial in the unknown T, a first-order newton method is used, applying the analytic gradient. Solution 
is facilitated by a small set of subroutines from GASP 16 which have been mex-ed for use with Simulink. 


Table 1. Coefficients used to compute saturation temperature from pressure 


Substance 

ii/io 1 

3 2/10 2 

h 

i4/io- 2 

is/io- 4 

ie/10- 6 

h/ io- 8 9 

Hydrogen 

-0.1059 

-0.3525 

0.3787 

-2.360 

8.553 

-16.18 

125.0 

Oxygen 

-0.5150 

-2.563 

0.2628 

-0.3622 

0.2652 

-0.1000 

0.1542 


For cryogenic He, supply pressure is computed from temperature and density. Figure 2 on the next pagea 
compares the GASP pressure values at 20K with data from the National Institute of Standards (NIST) 
chemistry web-book. 17 Figure 2 on the following pageb shows the time history of cryogenic He pressure 
during boost and upper stage firing. 


II.D.l. Volume and wetted area 


A sphere of radius r filled to a height h contains a volume equal to 


T/ ,2 ^ 

V = 7r rh - 7 r — . 


(55) 


The fluid covers a ‘wetted area.’ equal to A = 2irrh. Solving for fluid height as a function of volume is 
algebraically cumbersome, so instead, we solve Eqn 55 numerically using a Newton method. 
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(a) GASP pressure from density 


(b) Pressure time history 


Figure 2. Cryogenic Helium Sample Results 


III. Simulink Implementation 

The equations above were implemented in a Simulink model. Due to its size the overall block diagram is 
shown in the Appendix. The system decomposition in Simulink is by function rather than form. An attempt 
was made to decompose by form (physical components) but this resulted in algebraic loops. In this section 
we will tour the model block by block explaining similarities and differences. 

III. A. Tank Models 

The state derivatives in Eqns 23, 40, and 47 are numerically solved and integrated in the ‘LH2’ and ‘L02’ 
blocks which are expanded below. The Matlab function blocks ‘H2tank’ and ‘02tank’ provide solutions 
for the state derivatives, and include logic to choose between Eqn 23 (fuel explusion mode) and Eqn 40 
(recirculation mode). The current ‘H2tank’ implementation uses Eqns 23 and 40 to account for a mixture 
of H 2 vapor and He gas'in the ullage. 

III.B. Ullage regulators 

Ullage gas injection is modeled using a relay and the compressible flow relationship from Eqn 53. The relay 
is set to ±2.5 psi of the specified ullage pressure. The Matlab function block ‘GHel’ implements Eqn 53 
with He properties. The He regulator blocks use an external constant for supply pressure since the supply 
pressure for H 2 Helium injection is 435 psi (3.0 MPa) while the supply pressure for 0 2 He injection is 363 psi 
(2.5 MPa). The regulators for H 2 also include switching logic to select He injection for recirculation mode 
and gaseous H 2 injection for fuel expulsion mode. 

III.C. Relief Vents 

The H 2 tank vent relay is set to open at 49 psi (337.9 kPa) and close at 40 psi (275.8 kPa). Flow through 
the vent is determined using Eqn 53 where the gas properties are found by a weighted average of H 2 vapor 
and He gas properties according to the respective masses in the ullage. The 0 2 tank vent relay is set to 
open at 49 psi (337.9 kPa), and close at 30 psi (206.8 kPa). The flow through both vents is assumed to be 
sonic which is a valid assumption except in the case that the 0 2 tank pressure is below 30.1 psi (207.5 kPa), 
while the ambient pressure is standard sea-level or higher. 
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(a) Oxygen (b) Hydrogen 

Figure 3. Fuel and Oxidizer tank models 



Figure 4. Hydrogen ullage regulation block for fuel expulsion mode 


III.D. Recirculation Loop 


For recirculation modeling, the parts of the J-2X in contact with the fuel and oxidizer streams are modeled 
as a thermal capacitance. 


Tj2X = 


Q J2X 


m J2 XCj2X 

The heat transfer rate is modeled as forced convection. 


(56) 


Qj2x = hj2Aj2(Ti - T j2x ) (57) 

Fuel and oxidizer fluid warming are modeled using constant specific heat such that 

rp rp , Q J2X /rn\ 

T re -Ti-\ : (58) 

m re ci 

The simulink implementation for the O 2 and H 2 sides are identical except for fluid properties (including 
mass flow rates). Thus we show only the 0 2 half. Area Aj 2 , and' heat transfer coefficient hj 2 are external 
constants for trade studies and further refinement. 


III.E. Prevalve and Recirculation Isolation Valve 

The prevalves and recirculation isolation valves are simulated by simple time driven relays, and gains. Sensors 
for both open and closed positions are simulated by a boolean not. Flow rates are determined by a gain. 
Apart from flow rate gains, the 0 2 implementations are identical to those shown below for H 2 . 
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(a) Oxygen (b) Hydrogen 

Figure 5. Ullage regulation blocks for gaseous He injection 



(a) Oxygen (b) Hydrogen 

Figure 6. Relief Vent blocks for O 2 and H 2 tanks 


III.F. Liquid Level Floats 

The liquid level sensors are implemented directly from the liquid propellant / oxidizer mass remaining in 
the tank as shown below. Floats are simulated for overfill, 100% plus, 100%, 100% minus, 98%, 5%, and low 
level. The O 2 implementation is identical to that shown, with a different gain to reflect tank capacity. 

III.G. Tank Bottom Pressure 

Current vehicle acceleration is imported from a Maveric simulation of a nominal vehicle trajectory. Tank 
bottom pressure is determined by finding the liquid height in the tank, then multiplying by the fluid density, 
and current vehicle x axis acceleration. This value for fluid column head is then added to the ullage pressure. 
The liquid height in the O 2 tank is calculated by numerically solving Eqn 55 for h. The H 2 liquid height is 
calculated using a cylindrical tank model. 


T_J2-X 



Figure 7. Simulink implementation of the recirculation loop 
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(b) Recirculation Isolation 


(a) Pre valve 

Figure 8. Prevalve and Recirculation Isolation valve simulation blocks 


102 bottom press 




Figure 9. Tank bottom pressure sensor simulation blocks 


III.H. Cryogenic Helium Supply pressure 

The current value of cryogenic He pressure is computed by GASP from current density and H2 vapor 
temperature. The mass remaning in the cryo He tanks is found by integrating the net mass flow rates 
through the ullage regulators. The mass remaining is then divided by the total tank volume to find density. 
Here we show the simulink implementation. 



mdotHe02 

Figure 10. Simulink model of cryogenic He tank pressure 


IV. Simulation Results 

IV. A. Validation 

The basic function of our model was compared to an equivalent model in GFSSP 6 to show the validity of 
our model. Although the fundamental physics modeled in the two systems are essential the same,, there are 
some major differences in the numerical schemes involved. Most significantly, GFSSP uses an upwind finite 
differencing scheme such that the time marching equations are reduced to algebraic equations whose solution 
has accuracy similar to an Euler integration scheme. Both methods are currently using the GASP package 
for fluid properties which ensures a high degree of accuracy. Calculation of fluid to ullage and ullage to wall 
heat transfer coefficients were identical between the two methods. Other heat transfer coefficients are not 
well documented for GFSSP. 

Several factors influence the frequency of ullage repressurization oscillations perceived as the ‘sawtooth’ 
pattern of pressure response. The most significant at any time is the ullage volume. We adjusted the ullage 
volume in our simulation to roughly match the GFSSP response. Subsequent tuning was performed by 
choosing heat transfer coefficients from outside air to tank wall and tank wall to liquid propellant / oxidizer 
as adjustable parameters. These free variables were then tuned so that the simulink model pressure response 
matched the GFSSP response best in a least squares sense. The results are shown above. The time axes are 
zoomed above to aid in comparison. 

The validation results are promising. Although the oscillations eventually move out of phase with GFSSP, 
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(a) Hydrogen Tank (b) Oxygen Tank 

Figure 11. Validation using GFSSP 

note that the upward and downward slopes are quite comparable. The largest differences are caused by the 
GFSSP pressurization valve remaning open too long. Simulink automatically adjusts the simulation time 
step to precisely capture discrete changes such as valve openings / closings. GFSSP instead requires that 
the valve opens over a finite time interval, and does not check solution accuracy on the fly. 

IV.B. Detailed Results 

Figures 12 and 13 display simulation results for first stage boost (re-pressurization) and second stage firing 
(pressurization). Each pressure plot shows the partial pressure of ullage gas species, total ullage pressure 
and tank bottom pressure. Tank bottom pressure is essential to maintaining net positive suction head, and 
is restricted by structural limits. Each tank experiences a step in ullage pressure at the beginning of the 
simulation. First stage vehicle acceleration has an obvious effect on tank bottom pressure. 

The temperature plots show vapor temperature, liquid temperature and film temperature. It is important 
to remember that the film is modeled as an infinitessimal control volume, and thus film temp is assumed to 
be the saturation temperature at current ullage pressure. 

Ullage gas mass illustrates the key difference in Oxygen and Hydrogen tanks. Hydrogen vapor is injected 
into the Hydrogen tank to maintain ullage pressure during mainstage. Helium is injected into the Oxygen 
tank. 

Given that our results are comparable to those from GFSSP, we favor the use of Simulink in further 
development for the following reasons. 

1. Numerical Robustness: GFSSP’s numerical solution can have stability problems if the initial conditions 
(especially temperatures) are not set carefully. Simulink has built-in adjustable step algorithms that 
ensure accuracy of the solution. 

2. (Ex) portability: In GFSSP, the equations are hidden to the user. It serves as a good tool to create 
predictions, but would require extra effort to port the actual equations to another application say a 
fault detection algorithm 

3. Ease of adding features (and faults): GFSSP provides the flexibility to add new features like heat 
transfer to solids, but before long the user must write and compile FORTRAN subroutines to implement 
these features. In Simulink additional valves and logic are easily implemented by adding relays or 
interval tests. Faults may be introduced in' many cases by introducing erroneous limits on switches. 

GFSSP does account for the difference in temperature between added ullage gas and fuel / oxidizer vapor. 
We found the volume and rate of cooling of ullage gas to be the most influential factors in the frequency of 
ullage pressure oscillations. Adjusting these factors in GFSSP to match the simplifying assumptions in our 
simulation caused the quality of agreement shown above. 
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Figure 12. Simulation results for Liquid Hydrogen tank 
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Figure 13. Simulation results for Liquid Oxygen tank 
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V. Conclusions 


A simulink model of major components of the Ares I upper stage main propulsion system has been 
developed. Equations are based on conservation of energy for a small number of control volumes, and 
constraints from an isothermal ullage gas assumption, and invariant total tank volume. Thermophysical and 
transport properties of the working fluids are computed using the GASP 16 package which has been ‘mex-ed’ 
for use with Matalb / Simulink. Tank pressurization results during fuel / oxidizer expulsion are comparable 
to results from an established resource. This model, still under development has the potential to provide 
realistic nominal and fault simulation of sensors critical to an abort decision on the Ares I. 
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Tue, Nov 6, 2007 12:05 PM 


Subject: RE: Export control / ITAR question 
Date: Monday, November 5, 2007 11:28 AM 

From: Nelson, Michael A. (MSFC-ER21) <Michael.A.Nelson@nasa.gov> 
To: "Burchett, Bradley T. (MSFC)[]" <bradley.t.burchett@nasa.gov> 
Conversation: Export control / ITAR question 

Hey Brad, 

Absolute numbers are not allowed to go out on the ARES. The equations 
are generally available from already available public documents. So it 
is fine to show all the equations and relationships as you have done. 

It is just the absolute numbers about ARES that can not be made publicly 
available. 

So take the absolute numbers out and try to maybe show "example" ranges 

of numbers to get your point across without giving away the absolute 

numbers. You could zero out the bottom of a range and instead o f going 

from 6 to 40 it goes from 0 to 34 or something of that nature, for you | 

poul'd^muitipiy ( everything, ., : byy some/ ~,u npgb, 1 is he;d amount . and -make ^ a ^ngt^^-tPy 

rthe ac tual _ARE£f;; number-S--bu-t^Kavei^ r 

<be"e# :; mQdTfl6.dr , t o a void release _of._ARE.S_ ,_de_s ign de tai ls to the . public p er 5 j 

I T AR~~r'e's tr±c trpn s . 
v 

Hope this helps. 

Michael 

Original Message 

From: Burchett, Bradley T. (MSFC)[] 

Sent: Monday, November 05, 2007 11:16 AM 
To: Nelson, Michael A. (MSFC-ER21) 

Subject: Export control / ITAR question 

Mike, 

I am hoping to present the attached paper at the 2008 AIAA MST 
conference . 

Mike Watson asked that I confer with you regarding how ITAR might apply 
to the attached paper. The paper currently includes some numbers such 
as ullage pressures specific to Ares I. Where do we draw the line? 

Would it be okay as is, or should I remove the numbers specific to Ares? 
Are there other design details that need to be deleted as well 
(references to J2X , etc)? 

Thanks for any help you can provide, 

Brad 

Bradley T. Burchett, Ph.D., P.E. 

Associate Professor 
bradley . t . burchettOnasa . gov 

NASA VREO fellow, MSFC EV-43, Jun - Dec 2007 

"Force necessarily implies a mind. Thought moves mind and mind moves 
things. Force is the product of mind. None can, therefore, assert the 
reign of 'law' or 'force' without demonstrating the existence of the 
personal God of the Bible, which is the fact he seeks to deny... A 
process, necessarily implies contrivance, marks of design, and an 
intelligent Designer, which is God. Since a process cannot be the 
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